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Abstract 



O ' We propose a measurement of leading neutrons spectra at LHC in order to extract in- 

^ ■ elusive TT+p and 7r+7r+ cross-sections with high pt jets production. The cross-sections for 

these processes are simulated with the use of parton distributions in hadrons. In this 
work we estimate the possibility to extract parton distributions in the pion from the data 
on these cross-sections and also search for signatures of fundamental differences in the 



1^ ' pion and proton structure. 
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1 Introduction 

In recent papers P-[3] we have considered the possibihty of useing the LHC as pion- 
proton and pion-pion colhder. Here we continue to study the prospect of making unique 
measurements to extract cross-sections for vrp and vrvr interactions at TeV energies. 

Motivations for the present analysis are quite obvious. As one of the simplest QCD 
bound states and as the Goldstone boson of chiral symmetry breaking, the pion is a very 
interesting theoretical object: its structure carries important implications for the QCD 
confinement mechanism and the realization of symmetries like isospin in nature. It is also 
of practical importance for the hadronic input to the photon structure at low scales. The 
latter is connected via Vector Meson Dominance to the meson structure, which is poorly 
known and thus often replaced by the pion structure. 

The parton distributions (PDFs and GPDs) of the nucleons are now well determined 
by global analyses of the precise data for deep inelastic lepton-nucleon scattering, Drell- 
Yan and prompt-photon production. The recent summaries can be found in [l]-[7], in this 
work we use distributions from [8] integrated into PYTHIA [9]. The covered region is 

10"^ < X < 1, ~ 1 GeV^ < g^ < 10^ GeVl 

Unfortunately, determinations of the pion structure have made little progress over the last 
decade (see ^ or [ID]- [12] and refs. therein). They are based on old Drell-Yan and prompt 
photon data at fixed target energies and large values of the partonic momentum fraction 
X > 0.2. Many details are still based on pure theoretical assumptions, partially the precise 
knowledge of the nucleon distributions, and the use of different sum rules [12] that relate 
nucleon and pion distributions. In order to improve the situation, it has been proposed 
to measure the (virtual) pion structure at low x (down to x ~ 10~^) in deep inelastic 
scattering (DIS) and photoproduction with leading neutrons at HERA [13], [II]- Since 
the pion is by far the lightest hadron, its exchange dominates the p — ?■ n transition and 
it will almost be on its mass shell, particularly at small values of the squared momentum 
transfer t between the proton and the neutron. Analysis and references to the DIS dijets 
data are presented in [I3|-[IS]- Now we have pion distributions in the region 

10-^ < X < 1, 5 GeV^ <Q^ < 1.31 • 10^ GeV^ 

The parton model of QCD gives us the simple representation of a proton as a three 
quark state and the pion as a quark-antiquark one. This difference in the internal structure 
of nucleons and pions can be investigated experimentally. A long time ago it was proposed 
to measure forward-backward asymmetry and jet multiplicities in pp and np reactions with 
high pt events [IT] . The assymmetry can serve a clear signal that partons have different 
momentum distributions in protons and pions. It is possible to perform such analysis for 
TTTT reaction also. 

In this note we consider the production of leading neutrons plus inclusive dijet state, 
i.e. processes of the type 

pp^n jjX (1) 

and 

pp -^ n jjX n. (2) 



These processes may allow us to extract parton distributions in the pion in an unprece- 
dentally wide kineniatical region: 



10"^ < a; < 1, -10 GeV^ < Q^ < 10* GeV^ 



for yi up to 14 TeV. 



2 Extraction of tt p and n n cross-sections from single 
and double pion exchange measurements. 

In this section we give an outline of calculations of pion exchange processes with leading 
neutron production. Diagrams for Single (SvrE) and Double (DvrE) pion exchange pro- 
cesses are presented in Fig. [1] Factors F^, can be normalized to the low energy data [HI [19] 
and expressed as 



FA^.t) 



TT^pn 



-t 



2ht^l-2an{t) 



(3) 



167r2 {t-mlf 

where —t ~ (g ^ + m'^^'^)/{l — ^), the pion trajectory is aT,{t) = a^(t — m^) with the slope 
a'^ ~ 0.9 GeV~^. ^ = 1 — xl, where xl is the fraction of the initial proton's longitudinal 
momentum carried by the neutron, and G^o„„/{4:n) = G^+ /{Sn) = 13.75 [IS]. From 

recent data [19], [20], we expect b ~ 0.3 GeV~^. 
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Figure 1: Diagrams of Single (SvrE) and Double (DvrE) pion exchanges. 

Absorbtive corrections S, S2 are estimated in our model for high energy diffractive 
scattering. Details of calculations can be found in [Tl[2]. Final formulas for cross-sections 
look as follows 



Here 
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(T^p(s; {m.p,mj) ~ ff^_,p(S; {rrip, t}), a^^(s; {m^}) ~ a^,„,,^_t(^; {^1,2}), 



since the main contribution conies from pions with very low virtuahties \ti\ < 0.3 GeV^. 
We are interested in the kinematical range 0.01 GeV^ < |tj| < 0.5 GeV^, C,i < 0.4, where 
formulaes @ , (^ dominate according to [2T1 [22] . 

The next question is how to extract vr p and vr it cross-sections from the data on SttE 
and DttE? The exact procedure is similar to the Goebel [23] and Chew-Low [21] method: 
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The behavior of corresponding functions is shown in the Fig. |21 When t is equal to 
mass of the pion squared, we have no absorbtion at all {S = 1) and extracted cross-sections 
are independent on the model for rescattering corrections. 

The real situation is more complicated, especially from the experimental point of view. 
It is rather difficult to measure transverse momentum of a fast leading neutron, we can 
only get some restrictions on t from the acceptance of detectors. We propose to use the 
model dependent integrated method presented by formulas 
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Figure 2: Function from the expression (jS]) at fixed ^ = 0.05. The boundary of the 
physical region t ~ — ?7ip^^/(l — C,) is represented by vertical dashed line. 
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Figure 3: Rescattering corrections integrated with formfactors for y^ 
and a/s = 7 TeV (dashed): a) S'(s, ^); b) S^2(s, ^o)- 



0.9 TeV (sohd) 



Functions S and 5*2 are shown in the Fig. [31 Models for rescattering give us theoretical 
errors. If we have the data on p p and anti-p p total and elastic cross-sections, these 
uncertainties could be reduced to the errors of the data. For example, without LHC 
measurements at 10 TeV theoretical uncertainties can be estimated only from model 
predictions and can reach 20% for the most popular models. These errors are low for 
energies less than 1.9 TeV, since we have precise measurements from Tevatron. 

Our method (|9]) which used a very narrow t interval was applied to the extraction of 
71+ p total cross-sections at low energies [T] . It was shown in [T] that extracted points are 
close to the real data and four different model predictions which is a clear signal of the 
validity of our method. 



3 Pion exchanges with dijet production 

Since it is possible to extract vr p and vr vr cross-sections (total, elastic, Drell-Yan, direct 
photon or inclusive dijet production and so on) from the LHC data it should be possible 
to use such results to look inside the pion as we usually do it with proton and anti-proton. 
In this article we consider only the case of the inclusive dijet production as an example. 
Let us consider vr p or tt vr scattering with dijet production as a general process of 
the type /ii + /i2 — > jet jet X (see Fig. IHfor definitions). Momenta of particles can be 
represented as usual (for any -p = {po,P3', p})'- 
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In the case of vr p scattering we have 

Phi = P-K 

and for vr vr scattering 
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Figure 4: Process hi + h2 -^ jet jet X. 
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where 6* is the scattering angle in the CM frame of partons a and b. 

The basic formula for inclusive two parton production in the collinear approximation 
looks as follows 

dcrh^+h2-^cdx{M^f^) s sr^ ddab^cd,^ . 
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dz 



where fi/h{x) is the number density of parton i (quark, anti-quark or gluon) with the lon- 
gitudinal momentum x phfl in the hadron h. Renormalization and factorization scales are 
hidden in G and a. We can reconstruct momenta of final partons from jets measurements 
and then use our method (l9l)-(TTT]) to obtain combinations of PDFs in the pion and a pro- 
ton. Since we know proton PDFs from other experiments, we can extract combinations 
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Table 1: Leading order parton-parton cross-sections daab-^cd/dz = ^^Aab-^cd{s,i,u) 
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and d&ab^cd/dz can be found in the Table [TJ 

4 Experimental possibilities. 

We propose to measure jet production in single and double pion exchange reactions (([T]), 
([2])) at the LHC using the CMS detector [25]. Diagrams of the reactions ([1]) and ([2]) are 
shown on the Fig. [5] (a) and (b) correspondingly. 

In these processes two jets are produced in the hard n^p and tt+tt^ scattering which 
allows us to study parton distributions in pions in a still unexplored kinematical region. 
In this chapter we discuss perspectives of such measurements at 7 TeV. The Zero Degree 
Calorimeters (ZDCs) [26] can be used [H [2] to measure the leading neutrons. The ZDCs 
are placed on the both sides of CMS, 140m away from the interaction point. They have 
electromagnetic and hadronic sections designed to measure photons and neutrons in the 
pseudorapidity region \ri\ > 8.5. 

The Monte-Carlo generator MONCHER v. 1.0 [27] has been used for numerical simula- 
tion of the processes ([1]) and ([2]). This generator is developed by two of the authors of the 
article for SvrE and DvrE simulation specially. The kinematics of SttE and DvrE reactions 
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Figure 5: (a) and (b): diagrams for the processes of SttE and DvrE with 2-jet production. 
(c) and (d): diagrams for the processes oi pp inelastic interactions with 2 jets and leading 
neutrons production, which can imitate processes (a) and (b). 

are defined by the relative energy loss ^„ and the square of the transverse momentum t„ of 
the leading neutron. The vertex pvrj^^tn is generated according to the model described in 
the Ref. [Ill2]- PYTHIA 6.420 j9] is used for the vrj^^tP ~^ -^ generation for the single pion 
exchange and 'n'^irt'^^irt ~^ -^ generation for the double pion exchange. All background 
processes have been generated by PYTHIA 6.420. Diagrams of two background processes, 
pp inelastic interactions with 2 jets and leading neutrons production imitating signal, are 
presented on the Fig. O (c) and (d). 

PYTHIA 6.420 predicts 90.76 mb for the pp total cross section at 7 TeV. The inelastic 
part of this cross section, which is interesting for us as background, consists of 48.4 mb of 
minimum bias events, 13.7 mb of single diffractive events and 9.3 mb of double diffractive 
events. MONCHER 1.0 predicts 1.31^1.85 mb for SttE and 0.17^0.30 mb for DvrE 
processes at 7 TeV^J- Uncertainty in the cross sections comes from the different models for 
TT+p and TT'^TT'^ interactions. In this study we use the most pessimistic estimations from 
the Donnachie-Landshoff parametrization [28]. Therefore before any selections the ratio 
of the signal and background processes at 7 TeV looks as follows: 

CdttE : ctsttE : o"total = 1 : 7.6 : 530. 

In the rest of this paper we seek effective criteria to raise the signal/background ratio for 
single and double pion exchange. 

For SvrE we selected events with signal from neutrons in the forward or backward ZDC 
and with the absence of neutrons in the opposite one: 



N^ > & N^ = & ^n < 0.4 
N^ > & N^ = & ^^ < 0.4. 

For the DvrE we selected events with neutrons on both sides: 

K>^ & N^ > & ^n < 0.4 & C< 0-4. 



(22) 



(23) 



Here, N^ (iV^) is the number of neutrons hitting the forward (backward) ZDC, C,^ (^^) 
is the relative energy loss of the forward (backward) neutron. The signal to background 
ratio becomes equal to 0.22 for the SvrE and 0.76 for the DvrE after selections (|22l) and (!23l) 



""^ Cross sections for SttE and DttE are given for ^„ < 0.4. This is a kinematical bound of the model, 
see Ref. [1]. 
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Figure 6: (a) Multiplicity distribution of jets from the SttE (red(solid)) and pp inelastic 
events (black(dashed)) selected by fl22|) . (b) Multiplicity distribution of jets from the DttE 
(red(solid)), SttE (blue (dotted)) and pp inelastic events (black(dashed)) selected by 



respectively. For SttE all background events are produced in the inelastic pp interactions 
(minimum bias, single and double diffraction). For DvrE 20% of background are imitated 
by SttE and 80% come from pp inelastic interactions. 

FiglE] presents distributions of the multiplicity N^'^'^'^ of the jets from SvrE, DvrE and 
pp inelastic events selected by f l22l) (a) and by ( l23l) (b). Signal (red solid histograms) 
and background (blue dotted and black dashed) have rather different distributions of 
jSfjets ^ For the analysis we selected 2-jet events dominating in the SvrE and DvrE produc- 
tion(around 40%). 

Fig. [7] (a) and (b) shows pseudorapidity distributions in r]^'^^^ and transverse momen- 
tum pP-i.^ ^ of jets from the 2-jet SvrE (red(sohd)) and pp inelastic (black(dashed)) events 
selected by 0221) . The pl^ * distribution of the jets from signal (red (solid) histogram) shows 
a more gentle sloping behaviour comparing with jets from background (black (dashed) his- 
togram). It can be used for the further signal/background separation. Vertical lines on the 
plot (a) show the total Barrel, Endcap and HF acceptance of CMS. Plot (c) on the Fig. [7] 
presents the distribution of the sum of jets and neutron energies, Ejjn = Eji + Ej2 + En, 
for the signal and background 2-jet events selected by ( l22l) . The last right bin of the 
distribution, peaking at 7 TeV, corresponds to the exclusive production of jets. Events 
with Ejjn < 7 TeV come from the inclusive jets production. It is seen that 2-jet signal 
events are produced in the exclusive process dominantly (ratio of exclusive to inclusive 
production is approximately equal to 8.4). For the background, inversely, inclusive pro- 
duction of 2-jet events is more intensive (exlusive/inclusive ratio is approximately equal 
to 0.23). We also use this difference for the further signal/background separation. 

Fig. |8] shows the same distributions as Fig. [7] for DvrE (red (solid) histogram), SvrE 
(blue (dotted)) and pp inelastic (black (dashed)) 2-jet events selected by fj23|) . The dif- 
ference between signal and background in the pj-^ * distributions becomes more essential. 
Practically, there are no background events at pl^ * > 10 GeV, which can be used for the 
total separation of the 2-jet DvrE events from background. Almost all 2-jet DvrE events 
selected by ( p2|) are produced exclusively. The ratio of exclusive to inclusive production 
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Figure 7: Distributions in rj^^^^ (a), transverse momentum pl^ * (b) of jets and the sum of 
jets and neutron energies, i^jj^ = Eji + Ej2 + En, from the 2-jet SttE (red(sohd)) and pp 
inelastic events (black(dashed)) selected by (l22l) . The last right bin of the distribution (c), 
peaking at 7 TeV, corresponds to the exclusive production of jets. Events with Ejjn < 7 
TeV on the plot (c) come from the inclusive jets production. Vertical lines on the plot 
(a) shows Barrel, Endcap and HP acceptance of the CMS. 
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Figure 8: The same as on the FigjT] for the 2-jet DttE (red(solid)), SvrE (blue(dotted)) 
and pp inelastic events (black(dashed)) selected by (1231) . 

is equal ^14. 

Analysing distributions of jets for signal and background we suggest the following cuts 



for 2-jet SttE events selection and 



|^jcts| ^ 5 

pf' > 30 GeV, 



|^jets| ^ 5 

pf' > 10 GeV, 



(24) 



(25) 



for 2-jet DttE events selection, pj-^ '^-cut can be varied depending on trigger requirements 
and number of detected events to optimize signal/background ratio. 

Events of the reaction ([T]) selected by (!22l) fc (!24l) have ~15% of background from 2- 
jet pp inelastic production with leading neutrons. Events of the reaction ([2]) selected 
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by (!23l) fc (!25|) have ~3% of background from pp inelastic and SvrE production imitating 
signal. The additinal requirement 



6990 < J2 (Ejets + Eneutrons) < 7010 GcV, 



(26) 



selecting events with exclusive jets production, allows to suppress background for ([T| 
down to the level 6.5% and completely suppress the background for (j5]). 
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Figure 9: Distributions of events in (a) Invariant mass of (jj) system for SvrE events 
selected by (I22l) fc (!24l) fc(l26ll is shown by red (solid) line. 3% of background from pp 
inelastic events is shown by black histogram, (b) Invariant mass of (jj) system for DttE 
events selected bv (I23|)fc(|2ll)fc(l26|). 

Invariant mass of the 2-jet system produced exclusively in the reaction ([T]) is shown on 
the Fig. [9] (a) for the events selected by (I22 l) fcfl2^fc( l26|) . Invariant mass of the 2-jet system 
produced exclusively in the reaction (|2]) is shown on the Fig. |9] (b) for the events selected 
by (l23|) fc (!25!) fc (!26|) . Efficiency of the signal selection depends on p^^ '^-cut dominantly. 
With selections (!22|) fc( p^ fc( l26|) we save 2% of the single pion exchange events ([1]) and 
9% of the double pion exchange events (^ with (|23|)fc(|25|)fc(l26|). 



5 Discussions and conclusions 

We propose to measure reactions of single (SvrE) and double (DvrE) pion exchange with 
2-jet production at the LHC with CMS [25] using the ZDC calorimeter [26] to detect 
leading neutrons. Numerical simulation of reactions (II]) and (EJ) has been performed with 
MONCHER 1.0 [27] event generator. Background events from 2-jet pp inelastic interac- 
tions have been generated by PYTHIA 6.420 [9]. In this study we investigated effective 
criteria for selection of events in reactions ([1]) and ([2]) and estimated signal/background 
ratio. On the generator level of simulation the perspectives for such measurements look 
quite positive. 

For 2-jet DvrE events we can suppress completely the background from SvrE and pp 
inelastic interactions using a trigger for neutrons from ZDC (selections (123]) ) and proper- 
ties of the jets measured in Barrel, Endcap and HF of the CMS (selection ( 125]) ). The rest 
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of the 2-jet DvrE events after all selections, 9%, is equivalent to ~ 8 ^ 12yub (uncertainty 
is caused by different predictions for DvrE cross section). 

For 2-jet SvrE events we suggested selections (1221) fc fl21|) fc (l2Bl) which use only trigger 
requirements for neutrons from the ZDC and properties of jets. These selections suppress 
background from 2-jet inelastic events almost completely and save ~ 2% of the signal, 
which is equivalent to ~ 10 ^ 14/ib. 

The data accumulated by the CMS detector (more than 300 pb~^ at a time of the 
writing of this text) gives chances to extract millions of pure 2-jet SvrE and DvrE events, 
which are exclusive dominantly, for the detailed investigation of PDFs in the pion. 

From the theoretical point of view it would be very impressive if we had parton distri- 
butions in the pion in a still unexplored kinematical region, since pion is a fundamental 
"participant" of the strong interaction. Also comparison of PDFs in the pion and the pro- 
ton (anti-proton) can shed light on the mechanism of quark confinement and differences 
in the internal (quark, gluon) field structure of mesons and baryons. 
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